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Abstract — Many real-world systems such as manufacturing systems, transportation systems and logistics/distribution 
systems that play important roles in our modern society can be regarded as multicommodity flow networks whose 
arcs have independent, finite and multi-valued random capacities. Such a flow network is a multistate system with 
multistate components and its reliability for level (d;c), i.e., the probability that k different types of commodity can be 
transmitted from the source node to the sink node such that the demand level d = (d v d 2 ,...,d k ) is satisfied and the total 

transmission cost is less than or equal to c, can be evaluated in terms of minimal path vectors to level (d;c) (named (d;c)- 
MPs here). The main objective of this paper is to present an intuitive algorithm to generate all (d;c)-MPs of such a 
flow network for each level (d;c) in terms of minimal pathsets. Two examples are given to illustrate how all (d;c)-MPs are 
generated by our algorithm and then the reliability of one example is computed. 

Keywords — Reliability, limited-flow network, multicommodity, multistate system, (d;c)-MP. 

I. Introduction 

Reliability is an important performance indicator in the planning, designing, and operation of a real-world system. 
Traditionally, it is assumed that the system under study is represented by a probabilistic graph in a binary-state model, 
and the system operates successfully if there exists one or more paths from the source node s to the sink node t. In such a 
case, reliability is considered as a matter of connectivity only and so it does not seem to be reasonable as a model for some real- 
world systems. Many physical systems such as manufacturing systems, transportation systems, and logistics/distribution 
systems can be regarded as flow networks in which arcs have independent, finite, and integer-valued random capacities. To 
evaluate the system reliability of such a flow network, different approaches have been presented [7, 9, 14-23, 26-28].. 
However, these models have assumed that the flow along any arc consisted of a single commodity only. For such a 
flow network with multicommodity, it is very practical and desirable to compute its reliability for level (d;c), i.e., the 
probability that k different types of commodity can be transmitted from the source node to the sink node in the way that 
the demand level A = (d l ,d,,...,d k ) is satisfied and the total transmission cost is less than or equal to c. 

In general, reliability evaluation can be carried out in terms of minimal pathsets (MPs) in the binary state model case and 
(d;c)-MPs (i.e., minimal path vectors to level (d;c) [3], lower boundary points of level (d;c) [12], or upper critical 
connection vector to level (d;c) [7]) for each level (d;c) in the multistate model case. The multicommodity limited-flow 
network with budget constraints here can be treated as a multistate system of multistate components and so the need of 
an efficient algorithm to search for all of its (d;c)-MPs arises. The main purpose of this article is to present a simple 
algorithm to generate all (d;c)-MPs of such a network in terms of minimal pathsets. Two examples are given to illustrate how all 
(d;c)-MPs are generated and the reliability of one example is calculated by further applying the state-space decomposition method 

[4]. 

II. Basic assumptions 

Let G = ( N , A,U ) be a directed limited-flow network with the unique source s and the unique sink t, where N is the set of 
nodes, A = {a i \\<i <n] is the set of arcs, and U = (u v u 2 ,...,u,j ) , where U i denotes the maximum capacity of each 
arc for i = 1,2,..., n. Such a flow network is assumed to further satisfy the following assumptions: 

1. Each node is perfectly reliable. Otherwise, the network will be enlarged by treating each of such nodes as an arc 

[1]. 

2. The capacity of each arc tf is an integer-valued random variable that takes integer values from 0 to If 
according to a given distribution. 
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3. Every unit flow of commodity £ consumes a given amount p of the capacity associated with each arc. 

4. The capacities of different arcs are statistically independent. 

5. Flow in the network must be integer-valued and satisfy the so-called flow-conservation law [10]. This 
means that no flow will disappear or be created during the transmission. 

Assumption 4 is made just for convenience. If it fails in practice, the proposed algorithm to search for all (d;c)-MPs is 
still valid except that the reliability computation in terms of such (d;c)-MPs should take the joint probability 
distributions of all arc capacities into account. 

Since there are k different types of commodity within the network, the system demand level can be represented as a k-tuple 
vector d — {d x ,d d k ) where cl . is the demand level of commodity j for j = 1,2 ,...,k. Let X = (x { X 2 ,..., X n ) be a 

system-state vector (i.e., the current capacity of each arc CL under X is X ; , where X ( takes integer values 0,1,2,..., U j ), and 
V{X ) = (V ( X ) v V (X V(X ) l; ) , the system maximal flow vector under X where V( X) denotes the maximal flow 
of commodity j under X. (Whenever k > 2, there may be more than one maximal flow vector for each X. See the Appendix 
for more details.) Under the system-state vector X — (Xj X 2 ,..., X ) , the arc set A has the following three important subsets: 
N x = {a, g A I x,- > 0} , Z x = [a i e A I x ; = 0], and S x = {a t e N x \ V(X -e ; ) < V(X)}, where e t = (S n ,S j2 ,..., S in ), with 
S.. — 1 if j = l and 0 if j A i . In fact, A = S x kJ ( N x \ S x ) Z x is a disjoint union of A under X. 

A system-state vector X is said to be a (d;c)-MP if and only if: (1) its system capacity level is d (i.e., V(X)=d), (2) each 
nonzero -capacity arc under X is sensitive (i.e., N. = S x ), and (3) the total transmission cost is less than or equal to c. If 
level (d;c) is given, then the probability that k different types of commodity can be transmitted from the source node 
to the sink node in the way that the demand level d = (d l ,d 2 ,...,d k ) is satisfied and the total transmission cost is less than 
or equal to c, is taken as the system reliability. 

III. Model Building 

Suppose that P 1 , P 2 P'" are the collection of all MPs of the system, and let 

C = C k ,C 2 ,C 2 ,..., C 2 ,..., C 1 denote the transmission cost vector where C- is the unit transmission cost 

of commodity £ through arc CL . For each P 1 , W- = ^ { c' I a j e P J } and L. = min{» ( I <x g P 1 } are taken as the 

unit transmission cost of commodity i and maximum capacity through it, respectively. Under the flow- 
conservation law, any feasible flow pattern from s to t should satisfy that (1) the total flow-in and the total 
flow-out of each commodity for any given node (except for s and t) are equal, and (2) every unit flow of each 
commodity from s to t should travel through one of the MPs. Hence, under the system-state vector 

X = (Xj x,,..., X n ) with V(X) = (d l ,d 2 ,..., d k ) , any feasible flow pattern that the total transmission cost is less than or 

equal to c can be represented as a flow vector (//, f f*, f 2 , f 2 f* f f*) where fj is the flow of 

commodity £ transmitted through P J such that the following four conditions are satisfied: 


m 


= d , for each £ = 1,2,..., k 

i = i 

(1) 

k 

^ fjp £ < Lj for each j = 1,2,..., m 

£=l 

(2) 

k m 

for each i = 1,2,..., n 

£=\ 7=1 

(3) 

k m 

YL^fi 

t = i j = i 

(4) 
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Note that Hi/, : p l | e P' } is the least amount of capacity needed for a. under such a flow pattern 

t = 1 y=l 

k m 

(//./l 2 .-, fi’-’f!’ fm’-’ /m)» and so ’ under the system-state vector X, la,. eP j } 


/=! 7=1 


does not exceed the current capacity X : of tz ; . This fact is given in the following theorem. 


Theorem 1. Let X = (x v X 2 ,..., X n ) be any system-state vector for which V ( X ) = d. Then, the following is a necessary 
condition for the flow-conservation law to hold under X: 


k m 

{ fjp ' I a, e P '} for each i = 1,2,..., i 

7=1 j= 1 


(5) 


for any (/ t , /j ,..., /, , f 2 , f 2 ,..., f 2 ,..., f m ,f m f m ) which is a feasible flow pattern of flow d under X. 

Theorem 2. Let X be a (d;c)-MP. Then, the following is a necessary condition for the flow -conservation law to hold under X: 

k m 

for each / = 1,2,..., n (6) 

^=1 7=1 

for any (//, /, 2 ,..., f k , f 2 , f 2 ,..„ //,..., /J,,/j,..., /„*) which is a feasible flow pattern of flow d under X. 

fc m 

Proof . By Theorem 1, .r ^11 {/• p I a, e P 7 } for each i = 1,2,..., ?t. 

7=1 7=1 

1. For each e Z x , JC ( = 0 and so (6) holds. 

2. It remains to show that (6) holds for each a, e N x . Suppose, on the contrary, that there exists an arc <2 ; e N x such that 

k in ic m 

XE{/>' I a , (= P J } < Xf ■ Then, I fl, e P'} < jc f — 1 - In particular, V(X - e , ) - d = V(X) , and so 


7=1 7=1 


7=1 y=l 


A: m 


a t £ S x , which contradicts to the fact that X is a (d.c)-MP. Hence, x =XX( ia i * pi } for each a - e ^ • 

^=1 7=1 

The vector X = (x 1 ,X 2 ,..., X„) obtained by first solving F = (//, /j 2 ,..., //’ , f 2 ,f 2 ,..., f 2 / ( ‘, /^) subject 

to constraints (1) - (4) and then transforming such F = (// , f { 2 ,..., // , / 2 , f 2 ,..., f 2 ,..., /„', /„ 2 ,..., /j) to 
X = X 2 ,..., X n ) by applying the relationship in (6), will be taken as a (d;c)-MP candidate. To make it clearer that all 

(d;c)-MPs can be generated by the proposed method, the following theorem is necessary. 

Theorem 3. Every (d;c)-MP is a (d;c)-MP candidate. 

Proof. Let X = (x l ,x 2 ,..., X ) be any (d;c)-MP. By definition, we know that the maximal flow from s to t under X is d (i.e., 
V (X ) = d ) and the total transmission cost is less than or equal to c. Hence, under the system-state vector X, there exists at 
least one feasible flow pattern F = (f\ f 2 ,..., f k , f 2 , / 2 \..„ / 2 \..„ , /„ 2 ,..., f' m ) of flow d = (d v d 2 ,...,d k ) such that 

conditions (1) - (4) are satisfied. As X =(x 1 ,x 2 ,...,x n ) is a (d;c)-MP, we thus conclude, by Theorem 2, that 

k m 

x. = ^ y !./,'P' I £ P ' } for each i = 1,2,..., n . This means that X is a (d;c)-MP candidate. Hence, every (d;c)-MP is a 

7=1 7=1 

(d;c)-MP candidate. 
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In this article, we first find feasible solutions F = (f , / 



fm ) sub J ect t0 constraints 


(1) - (4) by applying an implicit enumeration method (e.g., backtracking or branch-and-bound [11]) and then transform such 
integer-valued solutions into (d;c)-MP candidates ( X, , X 2 , . . . , X n ) via the relationship in (6). Each (d;c)-MP candidate X 

must be checked whether all nonzero -capacity arcs under X (i.e., arcs N x ) belong to S x . If the answer is “yes”, then X 
is a (d;c)-MP. Otherwise, X is not a (d;c)-MP. The following two theorems play the crucial roles in checking whether a 
(d;c)-MP candidate is a (d;c)-MP. 

Theorem 4. For each (d;c)-MP candidate X, there exists at least one (d;c)-MP Y such that Y < X . In particular, X is not 
a (d;c)-MP if such a Y satisfies Y < X (where Y < X if and only if y t < X, for i=l, 2, n and Y < X if and 

only if Y < X and y i < X ; for at least one i). 

Proof. IfX is a (d;c)-MP, then Y must be taken as X. Suppose that X is not a (d;c)-MP; then there exists a nonzero-capacity 
arc a i (i.e.,a i s N x ) such thatV(X — e ; ) = V(X) = d . Let X 1 = X —e t ■ ■ Suppose that X 1 is a (d;c)-MP; then Y is taken 

as X 1 . Otherwise, the same procedure may be repeated for X 1 . However, this procedure will stop in finite steps, i.e., there 
exists an integer p such that X p < X p ~ l < ... < X 1 < X with V (X p ) = d and N xP = S x „ ■ The proof is thus concluded 

by letting Y = X P . 

Theorem 5. If the network is acyclic (i.e., contains no directed cycle), then each (d;c)-MP candidate is a (d;c)-MP. 

Proof. Let X = (x p X 2 ,..., X n ) be any (d;c)-MP candidate. By Theorem 4, we know that there exists a (d;c)-MP 
Y =(y l ,y 2 ,..., y n ) such that Y < X . Since V {X — Y) = d — d = 0, no flow is transmitted from s to t under 
X - Y = (x, — y, , x 2 — y 2 ,—, x n — y n ) . Hence, in case X ^ Y , I ={i\x j -y i > 0} is not empty and so [a { I i e 1} , which 
is a subset of N x , must form cycles since the flow conserves at each node (except for s and t) and there is no other sink 
except t (see Ford and Fulkerson [10] or Ahuja et al. [2] for more details). This means that if the network is acyclic, then 
I = (j) and so X = Y , i.e., each (d;c)-MP candidate X is a (d;c)-MP . 

Suppose that X 1 , X 2 ,..., X q are total (d;c)-MP candidates. We can thus conclude, by Theorem 4, that X J is a (d;c)-MP if 


X ' <j: X' for all j = 1 , 2 ,...., q but / ^ i . 

IV. Algorithm 


Suppose that all MPs, P l , P P"‘ , have been stipulated in advance [5-6, 24-25], the family of all (d;c)-MPs can then be 
derived by the following steps: 


Step 1. For each P J (j — 1,2,..., m) , calculate 



Step 2. Find all feasible solutions F = (ff / | 2 ,..., f \ k , ff f 



f 2 ?**•? fm ’ fm fm ) SUb J eCt t0 the following 


constraints by applying an implicit enumeration method: 


m 


(!) £//=^ for each £ = 1,2 ,...,& 


M 

k 


(2) y//y <Lj for each j = 1 , 2 ,..., m 


k m 



k m 
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where fj is a nonnegative integer for j = 1 ,2,..., m and l = 1,2,..., k. 

Step 3. Transform the solutions (//, fj, fj, / 2 2 ,..., / 2 \.„, /„', fj,..., /„*) into (d;c)-MP candidates 

fc m 

X = (x I ,x 2 ,...,x n ) via X ; la ; eP J } for i = l,2,...,n. 

t = 1 y=l 

Step 4. Check each candidate X one at a time whether it is a (d;c)-MP: 

(A) If the network is acyclic, then each candidate is a (d;c)-MP. 

(B) If the network is cyclic, and suppose {X 1 , X 2 X q } is the family of all such (d;c)-MP candidates, then X' 
is a (d;c)-MP if X J <j: X' for all j = 1,2,..., q but / ^ / . 


V. Examples 

The following two examples are used to illustrate the proposed algorithm: 


Example 1. 



Fig. 1. A SERIES -PARALLEL NETWORK. 

Consider the network in Fig. 1. It is known that U — (u l ,u 2 ,u i ,u 4 ,u 5 ) = (1,2, 2, 2,1) , C = (2, 3, 2, 3, 5, 6, 2, 3, 2, 3) , 
p = (yOj,/? 2 ) = (l,2), and there exists three MPs; P 1 ={a p a 2 },P 2 ={a 3 },/ )3 ={a 4 ,a 5 }. Given d = (1,1) and C = 12,the 
family of (d;c)-MPs is derived as follows: 

Step 1. L, = min{ l,2} = 1, L,- min{2} = 2, L 3 - min{2,l} = 1, W, 1 =2 + 2 = 4, W! 2 =3 + 3 = 6, W\ = 5, W 2 2 =6, 
^=2 + 2 = 4, 1E 3 2 =3 + 3 = 6. 

Step 2. Find all feasible solutions (fx , f{ , f 2 ■> f l ■> ) subject to the following constraints by applying an implicit 

enumeration method: 

{/;+/ 2 i +/3 l =i 

U 2 + / 2 2 + / 3 2 =i 

fl X 1 + /j 2 X 2 < 1 
</ 2 1 x1 + / 2 2 x2<2 
/ 3 1 x1 + / 3 2 x2<1 

/ 1 1 x1 + / 1 2 x2<1 
/ 1 1 x1 + / 1 2 x2<2 
</^xl + / 2 2 x2<2 
/ 3 1 x1 + / 3 2 x2<2 
/ 3 1 x1 + / 3 2 x2<1 

4// + 6/, 2 + 5/j 1 + 6 f 2 + 4/ 3 ‘ + 6 f 2 < 12 

where fj is a nonnegative integer for j = 1,2,3 and l = 1,2. 

Total feasible solutions are F l = (1,0, 0,1, 0,0) and F~ = (0,0, 0,1, 1,0). 
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Step 3. Transform such feasible solutions into (d;c)-MP candidates X = (x 1 ,x 2 ,x 3 , x 4 , X 5 ) via 

3 

Xj = ^ ^ { f \ p ' I a t g P'} for i = 1,2,. ..,5. Then X 1 = (1,1, 2,0,0) and X 2 = (0,0, 2, 1,1) are total (d;c)-MP 

i=\ j 

candidates. 

Step 4. The network is acyclic, and { X ^ , X ) is the family of all (d;c)-MP candidates. Since X' -jc X ' , 
X 1 = (1,1, 2,0,0) and X 2 = (0,0, 2,1,1) are total (d;c)-MP 

Example 2. 



FIG. 2. A BRIDGE NETWORK. 


Table 1 

Probability distributions of arc capacities in Example 2 


Arc 

Capacity 

Probability 

Arc 

Capacity 

Probability 

a i 

3 

0.60 

a 4 

1 

0.90 

2 

0.25 

0 

0.10 

1 

0.10 

a 5 

2 

0.80 

0 

0.05 

1 

0.15 

a 2 

2 

0.70 

0 

0.05 

1 

0.20 

a 6 

3 

0.65 

0 

0.10 

2 

0.20 

a 3 

1 

0.90 

1 

0.10 

0 

0.10 

0 

0.05 


Table 2 

Unit transmission cost on each arc in Example 2 


Arc 

Commodity 

Cost 

Arc 

Commodity 

Cost 

a \ 

1 

2 

a 4 

1 

1 

2 

2 

2 

1 

3 

3 

3 

2 

a 2 

1 

2 

a 5 

1 

2 

2 

3 

2 

3 

3 

3 

3 

3 

a 3 

1 

1 

a 6 

1 

2 

2 

1 

2 

2 

3 

2 

3 

3 


Consider the network in Fig. 2. It is known that U - (u l ,u 2 ,u 3 ,u 4 ,u 5 ,u 6 ) = (3, 2, 1,1, 2,3) , () = (p t , p n ,p } ) - (1,2,1), and there 
exists four MPs; P l ={a v a 2 },P 2 ={a v a 3 ,a 6 },P 3 = {a 2 , a 4 , a 5 ] , P 4 ={a 5 ,a 6 }. Given d = (1,1,1) and c = 16, the 
family of (d;c)-MPs is derived as follows: 
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Step 1. L, =min{3,2} = 2, L, = min{3,l,3} = 1, F, = min{2,l,2} = 1, L 4 = min{2,3) = 2, 

W/ =2+2 = 4, l¥, 2 = 2 + 3-5, W t 3 =3 + 3 = 6, W, =2 + l + 2 = 5, W; = 2 + l + 2 = 5, VP 2 3 =3 + 2 + 3 = 8, 

W, 1 = 2 + 1 + 2 = 5, W, 2 =3 + 1 + 3 = 7, W 3 3 =3+2 + 3 = S,W} = 2 + 2 = 4, W 4 2 =3 + 2 = 5, 
and W 4 3 =3 + 3 = 6. 

Step 2. Find all feasible solutions / 2 2 , / 2 3 , / 3 , / 3 2 , fhft’ ft) sub J ect t0 the following constraints by 

applying an implicit enumeration method: 

</, 2 +/ i 2 +/ 3 2 +/ 4 2 =i 


f\ X 1 + /j 2 X 2 + /j 3 X 1 < 2 
A* x 1 + / 2 2 x 2 + / 2 3 x 1 < 1 
' / 3 1 xl + / 3 2 x 2 + / 3 3 xl<l 
/ 4 1 xl + / 4 2 x 2 + / 4 3 xl<2 


// x 1 + / 2 * x 1 + /j 2 x 2 + / 2 2 x 2 + /j 3 x 1 + / 2 3 x 1 < 3 

// x 1 + / 3 x 1 + /j 2 x 2 + / 3 2 x 2 + /j 3 x 1 + / 3 3 x 1 < 2 

< / 2 x 1 + / 2 2 x 2 + / 2 3 x 1 < 1 
< / 3 1 xl + / 3 2 x 2 + / 3 3 xl<l 
/a x 1 + / 4 x 1 + f; x 2 + f; x 2 + / 3 3 x 1 + / 4 3 x 1 < 2 

fi x 1 + /i x 1 + / 2 2 x 2 + / 4 2 x 2 + / 2 3 x 1 + / 4 3 x 1 < 3 


4/; + 5/r + 6./I 3 + 5/2 + 5/ 2 2 + 8/ 2 3 + 5/3 1 + 7/ 3 2 + 8/ 3 3 + 4/4 + 5 f; + 6 / 4 3 < 1 6 


where /' is a nonnegative integer for / = 1,2, 3, 4 and £ = 1,2,3. 


Total feasible solutions are F 1 = (0,0, 1,1, 0,0, 0,0, 0,0, 1,0), F 2 = (0,1, 0,0, 0,0, 0,0, 0,1, 0,1), F 3 = (0,1, 0,1, 0,0, 0,0, 0,0, 0,1), 
and F 4 = (1,0, 1,0, 0,0, 0,0, 0,0, 1,0). 


Step 3. Transform such feasible solutions into (d;c)-MP candidates X =(x 1 ,x 2 ,x 3 ,x 4 ,x 5 ,x 6 ) 

x. = ££{/// \ a .e pJ } for i = 1,2,..., 6. Then X 1 = (2, 1,1, 0,2, 3), X 2 = (2, 2, 0,0, 2, 2), 

t= 1 i 

X = (3, 2, 1,0, 1,2) are total (d;c)-MP candidates. 


via 

and 


Step 4. The network is cyclic, and {X l ,X 2 , X 3 } is the family of all (d;c)-MP candidates. Since X' < X ' , every fd;c)-MP 
candidate is a (d;c)-MP. The result is listed in Table 3. 

Table 3 


List of all (d;c)-MPs in Example 2 


(d;c)-MP candidate 

(d;c)-MP? 

A 1 =(2,1, 1,0, 2, 3) 

Yes 

X 2 =(2,2, 0,0, 2, 2) 

Yes 

X 3 = (3,2, 1,0, 1,2) 

Yes 
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VI. Reliability evaluation 

If Y\Y 2 ,..., Y" ,,d:c) are the collection of all (d;c)-MPs, then the system reliability for level (d;c) is defined as 

*(d;c>= Priori* I* Sty 1 }}. To compute it, several methods such as inclusion-exclusion [8, 12], disjoint subset 

[13], and state-space decomposition [4] are available. Here we apply the state-space decomposition method to Example 2 
and obtain that R (d . c) = Prju”’] 1 '' 1 {X I X > Y' } } = 0.53235 for demand level d = (1,1,1) and c = 16 . 

VII. Conclusion 

Given all MPs that are stipulated in advance, the proposed method can generate all (d;c)-MPs of a multicommodity 
limited-flow network under budget constraints for each level (d;c). The system reliability, i.e., the probability that 
k different types of commodity can be transmitted from the source node s to the sink node t in the way that the 
demand level d = ( d l , d n d k ) is satisfied and the total transmission cost is less than or equal to c, can then be 

computed in terms of these (d;c)-MPs. This algorithm can also apply to the limited-flow network with single 
commodity. Hence, earlier algorithm [18] is shown to be a special case of this new one. 
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